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To ^ or not to that is the question! 
(adapted from J. Horvath's Concluding Remarks, with apologies to W. Shakespeare) 

We derive the form of the magnetic stress tensor in a completely general, stationary 
magnetic medium, with an arbitrary magnetization field M{r) and free current 
density j{r). We start with the magnetic force density / acting on a matter 
element, modelled as a collection of microscopic magnetic dipoles in addition to the 
free currents. We show that there is a unique tensor T quadratic in the magnetic 
flux density B(r) and the magnetic field H{'r) = B — AnM whose divergence is 
V • T = /. In the limit M = 0, the well-known vacuum magnetic stress tensor 
is recovered. However, the general form of the tensor is asymmetric, leading to a 
divergent angular acceleration for matter elements of vanishing size. We argue that 
this is not inconsistent, because it occurs only if M and B are not parallel, in which 
case the macroscopic field does indeed exert a torque on each of the microscopic 
dipoles, so this state is only possible if there are material stresses which keep the 
dipoles aligned with each other and misaligned with the macroscopic field. We 
briefly discuss the consequences for the stability of strongly magnetized stars. 



1 Introduction 

The equiUbrium thermodynamics of magnetized matter is important to un- 
derstand the structure and evolution of compact astrophysical objects such 
as white dwarfs and neutron stars. Of particular relevance is the pressure, 
since it enters in an important way in the equations that determine the struc- 
ture of the star (Newtonian hydrostatic equilibrium equation or its relativistic 
Oppenheimer-Volkov generalization) . If the gas that conforms the star is made 
of electromagnetically active matter (i.e., ionized matter or magnetic dipoles) 
in the presence of macroscopic electromagnetic fields, the total pressure will 
have both a kinetic and an electromagnetic component. 
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There seems to be some confusion and disagreement about the nature of 
the latter (a nd als o about the isotropy of the total pressure). For instance, 
recent paper^3^ claim that, for a magnetization parallel to the magnetic 
field, the electromagnetic contribution to the pressure takes the form 

P\\ =0, P^ = -MB, (1) 
where || and _L refer to the direction of the magnetic field. The authors of 
referencJ2l claim that, with an electromagnetic pressure of the this form and 
some assumptions about the growth of the magnetic field from the surface of 
the star towards its center, magnetars with surface fields of the order of 10^^ 
Gauss would be unstable to collapse. 

However, the correct electromagnetic pressure in the case of no magneti- 
zation, as derived from the Maxwell stress tensor for fields in vacuum is 

r>vac rjvac fn\ 

which clearly does not correspond to the zero magnetization limit of expres- 
sions (Q). 

In this work we address this problem by studying the electromagnetic 
force on a piece of magnetized matter from a purely classical point of view 
and constructing thereby the corresponding magnetic stress tensor for the sys- 
tem. Our results imply that the electromagnetic contribution to the pressure, 
defined as the negative of the diagonal components of the stress tensor, in the 
case of parallel magnetization and magnetic field is given by 

P„=-^, P.^^-MB, (3) 

which has the correct vacuum limit JJJ and, naturally, disagrees with ([T|l. 

Other works where the problem of the electromagnetic pressure has been 
touched upon arJSEEI. 



2 Magnetic force density on matter 

Macroscopic matter interacts with a magnetic field through its electrical cur- 
rents (due to motion of free charges) and macroscopic magnetization (due to 
the alignment of microscopic magnetic dipoles, usually associated with quan- 
tized spins). 

The force density (force per unit volume) on the currents is well-known 
to be f^^" = i X B/c, where j(r) is the electric current density, and B{r) is 
the magnetic flux density. 
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The force on the dipoles can be found by considering the force on a single 
dipole of moment m, whose Cartesian components are Fi = nijS/iBj. Adding 
the dipole moments in a small volume, we obtain the net force density, = 
MjViBj, so the total force density can be written aJ^ 

f^^^ eyfeij Bk + Mj V,Bj , (4) 

where Cy fe are the components of the totally antisymmetric (Levi-Civita) ten- 
sor, and the Einstein convention of summation over repeated indices is under- 
stood here and in what follows. 

We emphasize that this force density is not the same as obtained from 
defining a "magnetization current" j™"^ = V x A/ and writing f^°-9 = 
jmag y, B j although they agree when integrated over a bounded material 
body. This is most clearly seen in the simple case of a uniformly magnetized 
medium in a uniform magnetic field. Our expression for y^'""^ gives a uniform 
force per unit volume, as physically expected from the uniform distribution 
of dipoles acted on by the magnetic field. The alternative expression, for 
jmag ^ gjygg ^ vanishing force in the interior of the body (zero "magnetization 
current" in the interior), but a compensating force on the surface. 



3 The magnetic stress tensor 



3.1 Derivation 

The component Ty of the stress tensor is generally defined as minus the flux 
in the j— direction of the i— component of the linear momentum. The force 
density / on a matter element is the time derivative of its momentum density, 
i.e., minus the divergence of the momentum flux, thus 



/ = V • T, or 
Using Maxwell's equations 

V • S = and 



fi — ^ jTij ■ 



^ ri 47r ^ 
V xH=—j, 
c 



(5) 



(6) 



together with the definition H = B — AttM, eq. Q can be manipulated to 
write / as the divergence of a possible stress tensor, with components 



An 
BiBj 
An 



H,B. - [H ■ B --B' 



M.B, + M ■ BS,. 



(7) 
(8) 
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3.2 Uniqueness? 



Of course, the form just obtained is not unique, since any tensor with vanishing 
divergence can be added to T without changing eq. (Q. Thus, we search for 
a divergence- free tensor T, quadratic in the magnetic fields {B and H), whose 
general form is 

Tij — aijkiHkHi -f PijkiHkBi + jijkiBkBi, (9) 

where aijki , fiijki , a-nd 'jijki are components of constant tensors. These tensors 
must be rotationally invariant (in order not to introduce preferred directions 
other than those defined by the magnetic fields), which constrains them to be 
linear combinations of SijS^i, SikSji, and SuSjk- Thus, we write 

= {aH^ + PH ■ B + -fB^)6^j +a' H,Hj + (3' H,Bj + 13" B,Hj B,Bj, (10) 

where a, (3, 7, a', /3' , f3" , and 7' are constant coefficients to be adjusted so 
as to make the tensor divergence-free for any choice of the vector fields _ff (r) 
and B{r), with the only restriction that V ■ B — 0, i.e., VjBj = 0. 

In particular, we may first choose B = and H = Xaib (i.e., = XaSkt), 
where a and b are any given two different Cartesian indices, which implies 
\7jTij = 2aXaSia = (no summation over a), so a = 0. Changing H to be any 
other vector field with \7 ■ H ~ but {H ■ W)H / gives a' = 0. Exchanging 
the roles of H and B in the arguments just given, we obtain 7 = 7' = 0. 
Finally, putting H = XaBb and B = XcSd with different combinations of values 
for a, b, c, d (respecting the constraint that c ^ d to satisfy \7 ■ B = 0) implies 
/3 = /3' = P". 

Thus, among tensors quadratic in H and B, the one we originally obtained 
is unique in giving the correct magnetic force density. 



3.3 Vacuum limit 

In the absence of magnetization, the stress tensor given by eq. (O takes the 
form 

BiBj B'^ 
An 8n ' 

which is the correct, symmetric, well-knowrP form in vacuum. 



T^, = ^-^S.j, (11) 



3.4 Asymmetry and magnetic torque 

For arbitrary, non-parallel H and B, the tensor is manifestly non-symmetric, 
Tij 7^ Tji. It is usually arguec^ that the stress tensor must be symmetric. 
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on the ground that a non-symmetric stress tensor would produce a torque on 
a small matter element which decreases less quickly with decreasing volume 
than the moment of inertia of the element, thus producing a divergent angular 
acceleration. We note that this argument applies to the total stress tensor, 
which could be composed of (in principle non-symmetric) pieces of different 
physical origin. 

In the present case, the asymmetry is present only if H and B or, equiva- 
lently, M and B are not collinear. This means that the microscopic dipoles are 
not aligned with the macroscopic field, which therefore indeed causes a torque 
on each dipole, which add up to a torque per unit volume f = MxB. As a con- 
sistency check, we calculate a component of the torque density in terms of the 
stress tensor (as in Ref. 1^ but note the opposite sign convention for the stress 
tensor), — —eijkTjk — —eijkHjBk /{4:Tr), thus f = —H x B/^Air) = M x B, 
completing the consistency check. 

Of course, if there is no counter-acting torque (for example due to micro- 
scopic interactions among neighboring dipoles), the dipoles will orient them- 
selves along B, in which case the stress tensor becomes symmetric and the 
torque disappears. If there are microscopic interactions keeping the dipoles 
locally aligned and at fixed positions with respect to each other, then a macro- 
scopic matter element may act as a rigid body with a finite angular acceler- 
ation, that can not be subdivided into infinitesimal pieces with a divergent 
angular acceleration, because their magnetic torques are cancelled by the local 
microscopic interactions. 

3.5 Magnetized fluid and anisotropic pressure 

In a stationary, magnetized fluid, there are no microscopic forces keeping the 
the magnetization misaligned with the magnetic field, therefore Mir), B{f), 
and H(f) are collinear at every point f. Choosing a local Cartesian basis with 
z— axis aligned in the same direction, the stress tensor becomes diagonal, with 
components T^^ = Tyy = -P±_ = MB - B'^/{Stt) and T^^ = -P\\ = B'^/{87r). 
Thus, the parallel pressure P|| < always, corresponding to the usual tension 
along field fines, while Pi > as usual, unless M > B/{8n). 

4 Conclusions 

On the basis of physical arguments, we have found a magnetic stress tensor 
in magnetized matter that, although manifestly asymmetric, has the correct 
vacuum limit and appears to be consistent with all conceptual tests that we 
have applied. For the special case of a magnetized fiuid, with magnetiza- 
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tion parallel to the macroscopic field, we show that the longitudinal pressure 
remains negative and, unless the magnetization is extremely strong, the per- 
pendicular pressure remains positive. Thus, a neutron star with an extremely 
strong magne tic fi eld will definitely not collapse to a prolate structure, as has 
been proposed^^^ but, if anything, to an oblate structure, assuming that the 
star's fluid pressure does not impede it. 
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